Thermal Casimir drag in fluctuating classical fields 
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A uniformly moving inclusion which locally suppresses the fluctuations of a classical thermally 
excited field is shown to experience a drag force which depends on the dynamics of the field. It 
is shown that in a number of cases the linear friction coefBcient is dominated by short distance 
fluctuations and takes a very simple form. Examples where this drag can occur are for stiff objects, 
such as proteins, nonspecifically bound to more flexible ones such as polymers and membranes. 



The Casimir force, both quantum and thermal (or 
pseudo), arises due to the imposition of boundary con- 
ditions on quantum or thermal fields [IJ, the influence 
of the field is manifested by the force induced between 
two or more particles or surfaces in the field. However 
the presence of the field can also be seen by looking at 
the force exerted on a single particle when it is not at 
rest, for example a quantum friction exists for a neutral 
atom moving parallel to a dielectric surface and for 
parallel surfaces in relative motion [3|. While the lat- 
ter frictional forces require the presence of a dielectric or 
conducting surface, it has also been proposed that a fric- 
tional Casimir force can be induced by the uniform mo- 
tion of a polarizable molecule in a volume of blackbody 
of radiation which is in equilibrium in the rest frame of 
a containing cavity Q . Frictional forces acting on impu- 
rities in superfiuids and Bose Einstein condensates have 
also been studied, at low speeds no frictional force is 
found at a mean field level below a critical velocity 
however the scattering of quantum and thermal fiuctua- 
tions can be shown to generate a friction at arbitrarily 
small velocities [1^. 

For inclusions linearly coupled to classical fields driven 
by thermal fluctuations (for example a protein impos- 
ing a local curvature in a lipid membrane or a colloid in 
a binary liquid which is preferentially wetted by one of 
the phases) we recently showed that the insertion expe- 
riences a drag force when it moves at constant velocity 
(sj. The force arises because the insertion causes a local 
polarization of the fluctuating fleld, for instance a local 
magnetization for point-like fleld moving through a fer- 




FIG. f : Effects of inclusion interactions (thicker line) visual- 
ized on a one dimensional polymer (thinner line) (a) height 
constraint where the height of the polymer is fixed at zero to 
pass through a pore (b) a non flexible inclusion bound to the 
polymer which suppresses the local curvature. 



romagnet or a local non-zero curvature for a curvature 
inducing protein in a membrane. This local polariza- 
tion is not symmetric about the particle when the parti- 
cle moves at constant velocity and an effective frictional 
force is induced which opposes the inclusion's motion. At 
low velocity the force is typically linear in the velocity, 
f — ~Xv where A deflnes an effective friction coeflticient 
and at large velocities the force decays as |/| ~ 1/v. The 
average value of the force is in fact independent of the 
fluctuations of the field, it is mean fleld like and induced 
by the local symmetry breaking of the fleld by the in- 
clusion. An electrodynamical analogy for this linear cou- 
pling case is the drag on charged particles which create 
an average polarization (the polaron) in the surrounding 
medium. However a neutral atom does not induce an 
average polarization, but the fluctuations of the fleld in- 
duced by the dipole interact with the fluctuations of the 
fleld in the surrounding medium and it is this that leads 
to a frictional force. The soft matter analogue of this 
is drag experienced by an insertion in a classical fleld, 
driven by thermal fluctuations, which does not break the 
local symmetry of the fleld. In physical terms examples 
would be colloids in binary mixtures which are not prefer- 
entially wetted or membrane inclusion which flatten out 
membrane fluctuations, due to their stiffness, rather than 
inducing a local mean curvature. We show that a drag 
force is present for a variety of free fields undergoing 
stochastic dissipative dynamics of a type often used to 
model dynamics in soft condensed matter systems. 

The model we consider is for a scalar fleld with Hamil- 
tonian 
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= i/dx 



(x)A</>(x) + AH{^{z)) (1) 



where A is a positive deflnite operator which is model 
dependent, for instance for A (x—y) = [— V^-|-m^](5(x— y) 
the fleld (/> can represent the local magnetization in the 
Gaussian (high-temperature) approximation for the field 
theory of a ferromagnet [9|. Similar models can be used 
to represent the height fluctuations of fiuid-fluid/fiuid-air 
interfaces and lipid membranes [l^l and these models 
include operators with higher order derivatives due to the 
contribution of bending energy. The second term AH 
represents the energy of interaction between the field and 
an insertion at the point z. In this paper we will consider 
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interactions that are quadratic in the field variable (f>i the 
interaction with the insertion therefore does not affect the 
mean value of the field cj) about the inclusion (which is 
zero) but however modifies the field fluctuations. In this 
paper we will consider three interactions types 

AH = (a) ^02(z) (6) ^[V</.(z)]2 ic)^[V' cj^{z)]^ (2) 

For interaction (a) in the limit where h ^ oo the inser- 
tion field coupling induces a point-like Dirichlet bound- 
ary condition where 4> is energetically constrained to be 
zero at the inclusion's position z, an example of how this 
can occur is shown in Fig (la), a region of a polymer is 
pinned at a fixed height 0, this could be due to an optical 
trap or due to a pore that the polymer passes through. 
In the case (b) the gradient of the field is suppressed 
and in the limit /i — > cx) point like Neumann boundary 
conditions are imposed. Physically this corresponds to a 
region where the surface energy per unit area is large and 
the membrane is thus locally fiattened, this type of inter- 
action also occurs in the nematic phase of liquid crystals 
where inclusions suppress or enhance fluctuations of the 
director fleld about its average orientation [llj . Shown is 
Fig (lb) is the case (c), where stiff membrane inclusions 
suppress the local membrane curvature. 

We will study the force felt by the insertion when it 
moves at a constant velocity in the field. The dynamics 
of the field will be of a generic stochastic dissipative type 
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where the thermal noise driving the field is white in time 
and has a spatial correlation function obeying the local 
fluctuation dissipation relation 



(?7(x,i)r,(y,t)) =2T5(t-t')fi(x,y), 



(4) 



where T is the temperature. Above we use the opera- 
tor notation AB{x,y) = J dx' A{x,x')B{x' ,y), and we 
consider dynamical operators R which are invariant with 
respect to spatial translation, i.e. i?(x, y) = i?(x — y). 

We outline the computation of the force acting on the 
particle for Dirichlet type insertion/field interactions of 
type Eq. ([2^). The force on the particle is given by 



(5) 



where z = vi. The imposition of effective boundary 
terms via a local interaction with the field has the ad- 
vantage that the expression for the local force can be 
unambiguously written down as above for any field con- 
figuration ,13]. 

The average value of the force may be written in 
terms of the equal time correlation function of the field 
(,/.(x,i)0(y,t)) =C(x,y,t) as 



(f) = V,C(z,z') 



(6) 



In terms of the comoving coordinates of the insertion the 
correlation function can be shown to obey 

-C(x', y', t) ~ (v v., + V VyOC(x', y', t) 

= -(i?Ax' + RAy)C{x', y', t) - hR^,C{x', y', t)S{x') 
-hRyC{x',y', t)6{y') + 2TR{x' - y'), (7) 

where AxC(x, y) ~ J dx' A{x,x')C{x' ,y) denotes opera- 
tion on the (first) argument x of the correlation function 
C. In the steady state regime the double Fourier trans- 
form of the correlation function obeys 



i?(p)A(p) + i?(q)A(q) - *v . (p + q) ^(p, q 



dq' ~ 



dq' ~ 



+ f^^iP) J (^C(q',q) + /.i?(q) J (^C(q',p) 

^2T{27rfR{p)d{p + q). (8) 
In this notation the average force can be written as 

dq 



where 



(f) 



^(q) 



-ih 



{2n) 



q^(q), 



(9) 



(10) 



and can be shown to obey 

dp R{p)A{q) + R{q)A{p) 



^(q) 



{27:Y i?(p)A(p) -t- i?(q)A(q) - ■ (p + q) 
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A(q) 



(11) 



We can solve Eq. ([TT|) in two limiting cases to give a 
general picture of how the frictional force behaves. 
(i) Weak interaction limit: In the case where h is small 
we can compute A via a perturbation expansion in h, 
writing A ^ Aq + hAi + 0{h'^) we find that 

Ao(q) = J— (12) 



A(q) 



and 
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-T 



dp 



{2ttY 

R{p)A{p) + 7?(q)A(q) 



A(q)A(p) i?(p)A(p) + i?(q)A(q) - • (p + q) 



To this lowest order the force is given by 
Th'^ f dpdq 



(13) 



2 J (27r)2d 

a(P, q)(P + q) V ■ (p + q) (^(p)A(p) + i?(q)A(q)) 
A(q)A(p) ( fi?(p)A(p) + i?(q)A(q)l ' + [v (p + q)]^ 



(14) 
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with a(p,q) = 1. The equivalent of Eq. ([T4)) can also 
be derived for the case of the Neumann type interac- 
tion ~ Eq. ([5]d) - and curvature suppressing interac- 
tions ~ Eq. - and the corresponding results have 
o(p, q) = (p • q)^ and a(p, q) = p^q'^ respectively. In the 
above expressions the Fourier integrals may diverge and 
they are naturally regularized by an ultra-violet cut-off 
A = 7r/a where a corresponds to a microscopic size below 
which the field cannot fluctuate. In principal the cut-off 
could also be the size of the inclusion, however for sim- 
plicity in this study we consider inclusions of the same 
size as the microscopic cut-off. Here it is clear that for 
small V the force is generically linear in v but for large 
V we have! /I ^ 1/v as in the case of linear couplings to 
the field [8|. Finally we see that the force is fluctuation 
induced as it is proportional to T and is a second order 
effect in the coupling h. 

The small velocity regime: When v = wc find that 



^o(q) = T 



T 



1 



A(q)l + /^/(S.^ 



(15) 



For small v we assume that A{q) has a correction term 
Ai{q) which is of order v: 



Aiq) ^ Ao{q) + A,{q) + 0{v^ 
from which we find that Ai is given by 



(16) 



^i(q) 



r dp fl(p)Ao(q)+-R(q)Ao(p) 
., (2^^)" [fl(p)A(p)+_R(q)A(q)]" 

'^h^-^- r^-^ W77^^ • (17) 



{2rry fl(p)A(p) + fl(q)A(q) 



For simplicity we consider the limit of a Dirichlet point 
boundary condition where h —> oo and we find the result: 



A = 



cT 



d J dq' 



A(q') 



dqq 
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A(q)A(p) [fl(p)A(p)+fl(q) A(q)J 



R(p) 



i?^(p)A(p)+fl(q)A(q) 

(18) 

for the friction coefficient, with c = 1 and a = 1. In the 
same limit h oo we find a result of the same form with 
c = d and a ~ 2 for the Neumann case and c = 1 and 
a — A for the curvature suppressing case. 

As the operators A and R are isotropic, their Fourier 
transforms are functions of 9 = |q|, from Eq. (ITS)) and 
the forms for a given after, we find the result that all 
of these results are related as a function of the spatial 
dimension d via 

\N{d) = \D{d + 2)d and \c{d) AD(d + 4), (19) 

the subscripts D N and C denoting Dirichlet, Neumann 
and Curvature suppressing. 

Model A ferromagnet- Here we consider a simple fluctu- 
ating field where A(q) = q^ + rn? and we take simple dif- 
fusive dynamics for the field, characterized by a diffusion 



constant Dq so that R{q) = Dq. For d > 3 for Dirich- 
let and d > 1 in the Neumann case and all D for the 
curvature suppressing case, the friction coefficient takes 
a simple form in the limit where ma — > 0, i.e. where the 
microscopic cut-off a is much smaller than the correlation 
length ^ = 1/to of the field 



Az3(rf) = -j^Qid) 



(20) 



where Q{d) can be expressed as integrals independent of 
a and m. For example (5(3) ~ 0.479 in the case of a 
sharp cut-off as used here. However Q{d) will generally 
depend precisely on how the theory is regularized at short 
distances. In the cases d — 1 and d = 2 we find that 
the friction coefficient takes the form Xd = Tf{nia)/Do 
and where f{x) = — 21n(a;)/7r and f{x) — ln(— ln(a;))/2 
respectively. 

Model B ferromagnet- Here we again take the Gaussian 
ferromagnet but consider the case of model B dynam- 
ics which conserves the total magnetization and where 
R{q) = Doa^q^, where Dq is again a diffusion constant 
and the dynamical operator is given the correct physical 
dimensions via its dependence on the cut-off length scale 
a. In this notation we find for the Dirichlet interaction 
that for d > 4 A^i takes the form of Eq. I^U^ (the depen- 
dence on the cut-off a cancels with the factor of a^ ap- 
pearing in the definition of R and ultra-violet divergences 
in the integrals in the expression for A). For d = 3 we 

find that A Tln{ma)/DQ, A - T/[Doima\n{ma))'^] 

for d = 2 and for d = 1 we have A ~ T/[Do{ma)'^]. 
Rigid inclusion in lipid bilayers- Here we consider the 
case of a rigid inclusion in a lipid bilayer which locally 
suppresses the membrane curvature. The interaction 
term is of the form Eq. (Ht), and using the Helfrich hamil- 
tonian 12[ gives A{q) — Kq'^ + aq^ , where k and a are 
respectively the membrane bending rigidity and surface 
tension. Membrane fluctuation dynamics is dominated 
by the surrounding fluid and we have R{q) — l/Ariq, 
where rj is the viscosity of the surrounding fluid I4]. 
Here the drag is dominated by short distance behavior 
and we find the result Ac Trja/n, in the limit where 
a <C ^ = k/ct, (the correlation length of the membrane 
fluctuations) . 

Eq. (1111) can also be solved numerically for any h and 
V. We performed this numerical computation for a model 
A ferromagnet in d = 1; with a Dirichlet inclusion/field 
interaction and setting ma = 0.1, The computation was 
done with 101 Fourier modes. The mean drag force as a 
function of the velocity is plotted in Fig. ^ for different 
values of h. We see that for small v (/) is linear in v 
(and for large h the coefficient agrees with our analytical 
result) and that for large v it decays as 1/v. It is also 
interesting to note that the regime in v where linear fric- 
tional forces arise is very large, and its range increases 
upon increasing the interaction strength h. 
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FIG. 2: Drag force vs. velocity for a model A ferromagnet 
in d = 1 (on a log- log scale), with different values of h. The 
solid line is the analytical result Xdv. 



We have seen that in a number of cases, in particular 
when the insertion coupling depends on higher deriva- 
tives of the field and when the spatial dimension is high, 
the friction is dominated by the short wave length fluc- 
tuations of the field. In these cases it can be shown from 
Eq. (IT5|) (and dimensional analysis) that 



A - TAV(A) 



(21) 



where t(A) = [i?(A)A(A)]^^ is the relaxation time of the 
shortest wave length. The scaling of all the short-distance 
dominated results above can be deduced from this for- 
mula. The frictional force in this case is thus increased 
as the relaxation time of the short wave length modes 
is increased. The dependence of the friction coefficient 
on the microscopic cut-off is rather subtle. In the low 
dimensional ferromagnets with model A and B dynamics 
the friction coefficient increases as the cut-off decreases, 
but for curvature suppressing proteins in a Helfrich model 
membrane driven by stochastic hydrodynamics, the drag 
increases as the cut-off size increases. 

In the study of '3| of the drag on a polarizable molecule 
in blackbody radiation, one can see that the classical con- 
tribution to this drag is alsoproportional to T, as found 
in our study. However in [41 , the drag is of first order 
in the molecule's polarizability, whereas the effect found 
here is second order in the parameter h characterizing the 
interaction strength between the field and the particle. 
The drag experienced by irripurities in one-dimensional 
Bose Einstein condensates [7| shows a similar behavior 
to our result, in the high temperature regime the force is 
proportional to the temperature and perturbatively the 
force is second order in the coupling h of the impurity 
to the condensate. Physically the mechanism for drag 
creation is quite different in the dissipative systems stud- 
ied here, in the models presented here the inclusion is 



repelled from regions where the field fluctuates. The re- 
gion behind a moving inclusion fluctuates less than that 
ahead as the passage of the inclusion has flattened out 
the fluctuations. The field behind the inclusion is thus 
calmer than that ahead, thus making it energetically fa- 
vorable for the inclusion to move backwards and thereby 
inducing the drag force. 



The experimental observation of this drag requires 
a system where the fluctuation induced drag is not 
swamped by other drags, such as hydrodynamic Stokes' 
drag. Our result for drag on a curvature suppressing 
membrane protein appears to be only 1% of the hydro- 
dynamic drag given by the Saffman Delbriick formula 
[iSj for small insertions of the order of the membrane 
thickness. This is because the former is dominated by 
the viscosity of the embedding fluid and the latter by the 
membrane viscosity. However the addition of glycerol to 
the external solution can increase its viscosity by a fac- 
tor of up to 50 and in this regime the fluctuation induced 
drag will be of the same order as the hydrodynamic one 
and perhaps the drag predicted here could be observable. 
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